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Abstract 

For mechanical Weber gravitational wave antennae, it is thought that 
gravity waves are weakly converted into acoustic vibrations. Acoustic 
vibrations in metals (such as Aluminum) are experimentally known to 
be attenuated by the creation of electron-hole pairs described via the 
electronic viscosity. These final state electronic excitations give rise to 
gravitational wave absorption cross sections which are considerably larger 
(by four orders of magnitude) than those in previous theories which have 
not explicitly considered electronic excitations. 

1 Introduction 



Weber gravitational wave antennae are designed to weakly convert gravitational 
waves into acoustic vibrations P The damping of acoustic vibrations ulti- 
mately heat the detector accounting for the energy lost by the incident gravita- 
tional disturbance. The total cross section for the absorption of the gravitational 
wave is (of course) very small due to the weak nature of the gravitational forces. 
If g denotes an incident graviton which is converted into a phonon 0, then the 
first stage of graviton absorption is thought to involve the conversion 

9 - 0- (1) 

In the second stage of absorption, the phonon decays. In metallic samples (such 
as Aluminum) at not too high a temperature, the dominant mode of decay of 
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the phonon is into an 0I|S]electron-holc pair 

(t>~^e-+e+. (2) 

The decay in Eq. Q is known to be true in virtue of experimental work on sound 
wave attenuation in metals. The final electronic excited states are described in 
terms of an effective "electronic viscosity". Taking Eqs.lQJ and Q together, 
one finds that the gravitational wave decays into electronic particle hole pairs 
via the process 

g^^^e-+e+ (3) 

wherein the acoustic phonon (j) enters into the absorption process as a virtual 
(resonant) intermediate state. Also possible is the direct conversion of the gravi- 
ton into an electron hole pair 

g^e-+e+. (4) 

At first glance, it might appear that the electron-graviton coupling would be 
much weaker than the phonon-graviton coupling since the mass of the vibrating 
nuclei is much larger than the mass of the electrons. However, as will be dis- 
cussed in detail below, the graviton couples into the pressure and not directly 
into the energy density. The spatial (pressure) part of the space-time stress ten- 
sor Tij = Pij directly couples into a gravitational wave while the time (energy 
density) part of the space-time stress tensor Too — e couples into the Newtonian 
part of gravity 3 . 

The electron contribution to an earthly (ordinary) condensed matter pres- 
sure tensor Py dominates the nuclei contribution. The absorption rate of gravi- 
tational disturbances are considerably enhanced by the large kinetic energy and 
high density of final electronic states. A similar view has been discussed in some 
previous work0|7|. Nevertheless, the conventional wisdom is that the nature 
of the acoustic wave decay rate has little to do with the total absorption rate of 
the gravitational wave, provided that the former is much larger than the latter 
|S, ^ . Our purpose is to exhibit the flaw in the above conventional argument by 
providing a counter example. The fact that the electrons dominate the coupling 
to both the phonons and the gravitons considerably enhances the total cross 
section in metallic antennae (such as Aluminum). 

In Sec. 2, the gravitational wave propagator will be discussed. It will be 
shown that the "self energy" part of the propagator (as a fourth rank tensor 
"stress to strain ratio") can be viewed as a dynamical elastic tensor Young's 
modulus for distorting space-time. In Sec. 3, the gravitational wave will be 
treated as a weak curvature strain on a flat space-time background. The dy- 
namical elastic moduli of metallic antennae will be explored. In Sec. 4, the ampli- 
tude for scattering a gravitational wave off a metallic antenna will be computed 
in the mass quadrupole approximation. For solid metals, the final total cross 
section will be described in terms of the dynamical Lame elastic coefficient fi. 
In Sec. 5, the microscopic quantum mechanical expression for the gravitational 
wave total cross section will be exhibited. The mathematical derivations rely 
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on the Kubo formulae (in condensed matter physics) for the transport coeffi- 
cients associated with phonon damping. In Sec. 6, the dispersion relations for 
the dynamical elastic coefficients will be discussed. The results will be applied 
to the subtracted dispersion relations for the dynamic Lame coefficient and the 
viscosity. The electronic viscosity enhancement of the elastic Lame coefficient 
will be explicitly exhibited. In Sec. 7, the Feynman diagrams for the absorption 
of a graviton will be explored and the enhancement of the absorption due to 
electronic coupling will be computed. In the concluding Sec. 8, we present a 
brief summary of our results as well as directions for future work. 



2 Gravitational Propagators 

In relativity one considers (in general) a background curved space-time metric 
for the local proper time interval 

— c^dr^ — gf^iydx'^dx'^ (5) 

which arises from a solution of the Einstein field equations 

Upon varying by 5T^i, the "stress" applied to space-time 

T^i, —> T^i, + 5Tf_i^ (7) 

one finds a response in the distortion (i.e. "strain" {l/2)5g^i,) of space-time by 
requiring a neighboring solution to Eq.Q; i.e. 

9^il^ ^ gtiu + 5g^,v (8) 

The response in the strain due to the application of a stress defines the gravi- 
tational wave propagator D^ux^ via 

5g^,{x)= J D^,xa{x,x')6T^^x')dn' (9) 

wherein dfl = ^/—gd'^x. 

Let us decompose the stress into a source ST^^^ external to the antenna and 
a stress induced in the antenna ST^^ by the strain {l/2)5gfj_^. In total 

ST^- ^ST,^^, + ST,i-^. (10) 

The full propagator describes the response in the strain to an "external" stress 

Sg^,,ix)= [ A^,x4x,x')ST^^^t{x')dn', (11) 
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while the "induced" stress responds to the strain via the elastic response function 
6T,'a^') = Jn'^"^ix',x")SgM^")dn". (12) 
We note in passing that EQS.(f5|)- (fT^ imply Dyson's Ea. (fT5|) below; i.e. 

' D^,^p{x,x')Il<'P^\x',x")A^nXa{x",x"')dD.'d^"{\Z) 



wherein the dynamic elastic coefficients 11"''^'' play the role of the gravitational 
wave propagator self energy. 

Since the external stress radiates the incoming strain via the relation 



Sg'i^K^) - j D^,U=^,x')5T^^^{x')dn', (14) 

it follows from Eas. (|ll|) . and (|14|l that the gravitational wave scattering 
equation takes the form 

5g^,{x) = 5g^::\x) + f f D^,^p{x,x')TrP^'' {x' ,x")5gxa{x")dQ.'dQ," . (15) 



In Eg. pSfl . the total gravitational wave is the sum of an incoming gravitational 
wave and an "elastic" scattered gravitational wave 

5g^..{x) = 5g^:\x) + 5g'^/^t\x) 

Sg'^^Hx) - [ [ D^,^pix,x')U''f^^x',x")SgUx")dn'dn". (16) 



From the imaginary part of the forward elastic scattering amplitude one may 
compute the total cross section for the wave to scatter off the antenna. Since 
the propagator D^^af3 contains the weak gravitational coupling strength G to 
first order, it is sufficient to employ the lowest Born amplitude for the scattered 
gravitational wave 

^^(e;,Bo™)(^) ^ / I D^,^p{x,x'WP>^%x\x")5g'^:\x")dn'dil". (17) 



To compute the scattered wave to lowest order in G one needs to know (i) the 
incoming gravitational wave Sgi^^\ (ii) the dynamic elastic coefficients 11"'^'^°' 
of the antenna and (iii) the unperturbed (by the antenna) gravitational wave 
propagator D^^ap- Let us consider this in more detail. 

3 Gravitational Wave Antennae 

Consider a gravitational wave moving through a flat space-time background in a 
reference frame in which an antenna is at rest. In the transverse traceless gauge 
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one may describe the gravitational wave by a spatial strain tensor Uij (r, t) which 
enters into the proper time as follows 

c^dr^ = c^dt^ - ds^ 

ds^ = |drp + 2dr • u • dr. (18) 

The definition of "strain" in terms of the spatial metric in ds^ conforms to the 
standard usage in elasticity theorvfTU]. In terms of the traceless part of the 
spatial pressure tensor pij (r, t) in the antenna, 

p = P - il(irP), (19) 

the gravitational wave equation (with a pressure source) reads 

The general solution of Eq. (|20|l may be written 

u,M) = u«(M)+(5)/£fe^i^.V (21) 

where i? = |r — r'|. The incoming gravitational wave obeys 

For a dynamical elastic antenna, the linear causal response in the pressure 
p to a material spatial strain u is described by a non-local (in space and time) 
expression of the form 

p,j{r,t) = - J g,jkiir,r',t)uki{r',t-t)dtd^r'. (23) 

If the strain is turned on at a complex frequency ^ with 3toC > and if 

u(r,0=5ie{e-<*u(r;C)}, (24) 

then the pressure response 

p(r,t)-3?e{e-^';*p(r;C)} (25) 

is described by a frequency dependent Young's modulus y. The explicit expres- 
sion for the pressure is 

K,(r;C) = - / y^M^,r'■,C)uklir';Od^v' (26) 
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where the dynamical Young's modulus obeys 

roc 

3^„fei(r,r';C) = / e'^'g,,M{r,r' ,t)dt. (27) 
Jo 

The scattering equation follows from Eas. (|21|) and (|24ll - 12()f) : It is 

"..(r;C) = ug"^(r;C)+u!f (r;C) 
ui'\r-X) = -^||(^^)3^,,.Kr',r";CKKr";C)rfVdV' (28) 
where R — |r — r'|. 

Note that Eq. H28|l is the flat space-time limit of Eq. H16|l in thinly disguised 
form. Lowest order perturbation in gravitational coupling constant G yields 

(eLBorn) _ f f f f^' ^'f . /■\„,{in) f,f . ..,3/ ,3^f' 



<--°™'(r;C)^-^y J (^^j3^.,M(r',r";C)<r(r";C)dVd^r" (29) 

as the flat space-time limit of Eq. 117|) . Far away from the antenna, the scattered 
gravitational wave for real frequency obeys 



where 



F,, (u;) - - (^^^ J J e-*^ --3^,,feKr,r';u; + zO+)47^(r';a;)d3rd3r' (31) 

with the outgoing wave vector = (ujr/cr) ~ {ljv/c). If the incoming wave 
has the form 

v!f(v-u;)^{e,i),e^-\ (32) 
then the elastic scattering amplitude 

{f\H^)\i) = {e.,)}F,,{u;) (33) 

follows from Eq.lEIl) and ^ to be 



{e-*/--(e,,)}X,fe,(r,r';c. + zO+)(efe0.e''^'-""'} • (34) 

The total cross section ai for a graviton (in state i) to be absorbed by the antenna 
is determined by the imaginary part of the forward scattering amplitude 

a.(c^) - (^^) ^mii\J^iij)\i). (35) 

Eas. (|34l) and (|35() constitute the central result of this section. The total cross 
section a is determined by the dynamical elastic Young's modulus y of the 
antenna. 
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4 Mass Quadrupole Antennae 

If the wavelength of the incident graviton is large on the scale of the size L of 
the antenna, 

{io/c) « (1/L), (36) 

then one may ignore the wave vectors k^j in the integrals of Ea. (|34|l . This 
constitutes the mass quadrupole approximation. For an antenna of volume O 
and with 

^lx^MC,^)^ I I y^Jklir,r'■,C)d''rd^r', (37) 
In Jn 



the mass quadrupole absorption cross section is 

<J^{^) « - ^r^^{{e,,)*x^M^ + iO+; n){eki),} . (38) 

Averaging over the initial gravitational wave polarizations (i.e. helicity), 

and defining the transverse dynamical elastic Lame coefficient as 

^l{(,n;n) = jJ2 {i^^Jy^'X^MC:mekl)^'} (40) 
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yields the total cross section 



a(uj, n;n) = -( IMl^ ) <^mfi{LU + iO+ , n; n), (41) 

\ C-^iiJ J 

where n denotes a unit vector along the incident propagation direction of the 
gravitational wave. We note in passing that the dynamical Lame coefficient 
/x(C, n; f2) for a finite volume system may have frequency ^-poles corresponding 
to both transverse and longitudinal acoustic modes. The reason is that when 
a transverse sound wave reflects off the boundary surface, a longitudinal sound 
wave may be produced with finite amplitude and vice-versa. 

We now presume that a large antenna can be represented as a poly crystalline 
"isotropic" elastic system. The absorption cross section shall nevertheless de- 
pend on the incident propagation direction of the gravitational wave if the shape 
of the antenna is not spherical. 

The imaginary part of the Lame coefficient may also be expressed in terms 
of the real part of the viscosity coefficient according to 

- 9m {n{uj + iO+, n; n)} = uj^e {rj{uj + iO+,n; n)} . (42) 

Thus, with the mass of the antenna denoted by M = p^l, Ea. (|41|l reads 

a(., n; n) ^ ('-^) ^^M^ + . (43) 
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Ea. l|I5|l is the central result of this section. The absorption cross section for 
a graviton is directly related to the "kinematic viscosity" (ji/p) whose value 
determines the life time of the converted transverse phonon. 

5 Microscopic Expressions 

If one writes the microscopic dynamical Young's modulus as 

y^jki[vr'- = y.,,ki{^v'- 0) - icn,M-^v'; c), (44) 

then the Kubo formula for the non-local (in space) and frequency dependent 
viscosity response function is given by 

7i,M^v'-0 ^11^ [I e*'^* {pkl{r',-iX)p^,ir,t)) dtj dX, (45) 

where 



The viscosity is then 



From Eas. ((T^ . and it follows that the pressure integral 

Pd^r =\ j rrdivdivPd^r. (48) 

is of central importance. Ea. H48|) follows from twice integrating by parts. On 
the other hand, we have from conservation of mass 

^ + divJ = (49) 
and from conservation of momentum 

— + divP = (50) 
that ^ 

= divdivP. (51) 
From Eqs.gHI) and (|3TJ it follows that 

J P(r, t)d^r = i (^ly I rrp(r, t)d^r. (52) 
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For the problem at hand, the mass quadrupole of the antenna may be defined 

as 

D(<)= [ {3rr-r^l) p{r,t)£r. (53) 

From Eqs.lUni), lE2Jl and ^ it foUows that 

p(r,t)d3r=iD(i), (54) 
n D 

while Eqs.lQSl ^-i^d (ES imply 
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n Jfi 

[3 ( i-oc 



(DM{-i\)D,,{t))dt\d\. (55) 



iQhQ. Jq 

Finally, from Eas. l|17jl and (|^ it follows that 

v{C,n;n) = i ^(e,,)*,/i,,fe,(C,n;f})(e,0.'- (56) 



Eas. (|53|) . H55|l and H5t)|) yield a rigorously exact microscopic expression (to lowest 
order in G) for the total gravitational wave cross section in Ea.()43|l. 

6 Dispersion Relations and Sum Rules 



The viscosity coefficient obeys the analytic dispersion relation 
'2iC\ ^e{r]{Lj + iO+,n;n)}du; 



,{C,n;n) = - ( ^ ) / --ViV-^J^^-^-'f- for ^mC > 0. (57) 



With a finite viscosity ?7(0, n; il) and shear modulus /x(0, n; il), the full dynamical 
Lame coefficient reads 

m(C, n; n) = ^l{0, n; n) - tCviC, n; ft). (58) 
In the hmit |^| oo, Eas. (|S7|l and imply 

-/ 3fie{77(tj + iO+,n;f7)}dtj= (/i(oo,n;f7)-^(0,n;r2)). (59) 
Jo 

The microscopic evaluation of the sum rule follows from Eas. (|55|l . H56|) and 
H59|) via an equal time commutator 

Afi — ^(cx), n; il) — /^(O, n; fi) 

= — 5Re / (ey)*,/iyfci(w + iO+,n;0)(efci)i'rfw 

^ I fw) E (^'^■)*' ( [a, , ) (e.O.' I • (60) 



Ef^)+E-^^(-«^)' (62) 



a<b 



From Eq. I|53|) the mass quadrupole quantum operator may be written as 

D,,=J2Ma{3ra^ra,-6,,rl). (61) 

a 

wherein Ma is the mass and Va is the position of the a*^ particle in the gravita- 
tional antenna. These particles include both the electrons and the nuclei in the 
condensed matter. If we employ the Hamiltonian H = K + U which is the sum 
of the kinetic energy K and the potential energy U, i.e. 

Pi 

2M, 

a ^ 

then the time rate of change of the quadrupole operator in Ea. (|61() reads 

Dtj = ^ {i{raiPaj +Pairaj) " % (l"a ' Pa + Pa ' ^a)} ■ (63) 
a 

Similarly, 

-2^.E{(f)+(r-u}^ (64) 
where for the case of two particle potentials 

Pa = fa -Vat/ = "Vq ^ Uab{rab)- (65) 

b^^a 

The equal time commutation relations are computed employing 

[Pa%,rb:j] = -ihSabSij. (66) 
A very tedious computation employing Eas. HfciU|) - (|66|l yields the sum rule 
2 f°° 

Ap. = - ^efj{uj + iO+;n)dLLj 

Jo 

^ ik ^^(^iK)+'^(j2j-b<birab)^ + (j^J^-bKbirab)^^ 



(67) 



where f]{(^;fl) represents a spherical average over n in ri[(^,n;fl). For the case 
of a model Hamiltonian for which the potentials are completely described by 
Coulomb's law, 

U^e^Y.—' (68) 



a<b 
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Eqs.lEI) and ^ imply 

A/2 - - r dief]{u; + iO+: n)dcj = ( illMzlM^ . (69) 

From the virial theorem for a collection of particles (electrons and nuclei) inter- 
acting with pure Coulomb potentials, the thermodynamic pressure obeys 



while the energy per unit volume 



Eqs.(En|)-(|ZIll imply 



Afi=^J^ dief]{uj + iO+;n)du;^ (^^^^j^] . (72) 



Eqs.GOIl-GlJ) imply that 

A. ^(i|) <--(!). 

where the kinetic energy per unit volume 

t={K)/n. (74) 

For pressures of the order of one atmosphere (or less) one finds that t >> P. 
Thus 

wherein p and Vs represent, respectively, the mass density and the shear sound 
wave velocity. Note that the dynamic Lame coefficient /i(C, n; fi) contains C- 
poles at modes determined by both the longitudinal and transverse acoustic 
phonons. Nevertheless, the zero frequency static Lame coefficient p, = pVg 
depends only on the shear (transverse) wave sound velocity Vs- 

The mean kinetic energy per unit volume is dominated by the electrons. If 
n represents the number of electrons per unit volume, then a Thomas-Fermi 
uniform electron distribution estimate of the kinetic energy per unit volume of 
the electrons reads 



For example, an Aluminum atom has Z = 13 electrons and an atomic weight 
of A = 27 so that the number of electrons per unit volume is given by n = 
{Zp/AM). Here, M is the atomic mass unit. Thus, we find the estimate 



84 (Aluminum). (77) 
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Figure 1: Shown above are the foUowing processes: (i) A graviton decays into 
a particle hole pair g —^ e~ + e"*". (ii) A phonon decays into a particle hole 
pair (/) — > e~ + e"*". (iii) A graviton again decays into a particle hole pair, but 
there is a conversion to an intermediate resonant phonon g ^ cf) + . 

The conversion coupling strength is /2. (iv) A graviton decays into a particle 
hole pair, again with a conversion to an intermediate resonant phonon g 
[e~e~^] — !■ — > e~ + e"*" employing an electron loop coupling. The enhanced 
loop conversion coupling strength is A/i. 



The electronic contribution A/2 to the dynamic Lame coefficient leads to an 
electronic renormalization of the conversion of a graviton into a phonon as will 
be discussed directly below. 



7 Enhancement Factor for Graviton Absorption 

The Feynman diagrams for the physical processes here considered are shown in 
FigH The absorption rates are controlled by the material viscosity rj which for 
a metal (at reasonably low temperatures) is determined by the motions of the 
electrons, (i) If there is one graviton inside of a very large material object of 
volume rj, then the transition rate per unit time for a low frequency graviton 
to be absorbed is equal to the flux per unit time per unit area times the cross 
section; i.e. 

' C(7{uj) " 



r. = lim ( ^ ) (78) 



which, in virtue of Ea. H43|) . reads 

Fg = (^^^) V for g^e- + e+. (79) 

(ii) Similarly, the absorption of a phonon is also controlled by the material 
transverse and bulk viscosities, respectively, rj and r]b. For a phonon with wave 
vector k the transition per unit time for absorption by the material is given by 

r ch( longitudinal) — ) ?7b H ) for cj) ^ e + e^, 

\ P J \ 3 7 
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r ^{transverse) = ( j rj for — > e + e'^ . (80) 

(iii) In gravitational wave antennae the absorption of a graviton takes place 
through the intermediate conversion of a resonant phonon via 

g^0^e"+e+. (81) 

The matrix element for the conversion g — > is described by the Hamiltonian 

H'{g,<j>)^2fljis:u)£r, (82) 

where (in the antenna) s is the strain due to the phonon and u is the strain due to 
the graviton. Ea. H82|l describes the graviton conversions into both longitudinal 
and transverse phonons. Nevertheless, only the Lame coefficient fj, enters into 
the coupling. The reason is that the graviton is purely transverse and traceless 
no matter which kind of phonon it produces, (iv) An electronic enhancement 
of the graviton to resonant phonon conversion process arises from the electron 
loop process shown in Fig.l^; i.e. the one loop conversion process 

g -> [e-e+] ^ (83) 

determines the coupling strength renormalization 

fi^fi + Afl, (84) 

wherein A/2 has been evaluated in Ea. l(75|l . 

Since the absorption rate for gravitons is computed from the absolute value 
squared of the matrix element, the central point of this work is finally evident. 
The enhancement of the graviton to resonant phonon detection process via 
electronic loop viscosity is given by 



£ = 



(85) 



For the case of an Aluminum gravitational wave antenna, Eas. H77|l and H85|l 
imply a large electronic induced detection enhancement; It is 

£ « 7.2 X 10^ (Aluminum). (86) 

Such an enhancement allows for the experimental feasibility of detecting gravi- 
tational radiation from sources within our own galaxy. 

8 Conclusions 

In this work, we have considered in detail the dynamical response produced in 
a metallic Weber bar by a gravitational wave. Our treatment differs from all 
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previous analyses in several major aspects. First, departing from tradition, we 
follow the Einstein equations so that a (transverse, traceless) gravitational wave 
is directly coupled to the (transverse, traceless) part of the pressure. Secondly, 
observing that in a metal (at not too high a temperature) electrons (not the 
phonons) provide the bulk of the pressure. Previous estimates of the couplings 
based on direct graviton to resonant phonon conversions need a major revision. 
In fact, we find that the electron-hole pair loop contributions significantly renor- 
malize the relevant couplings. We have computed this using a Thomas-Fermi 
model with uniform electron density to estimate the average kinetic energy per 
unit volume of the electrons. The enhancement factor is found in this way to 
be £" ~ 7 X 10'^. The true factor is likely to be somewhat higher since the as- 
sumption of uniform electron density underestimates the kinetic energy of the 
localized core electrons. 

Our results can be employed to explain the solution to a theoretical problem 
raised by experimental observations with the room temperature resonant detec- 
tors of the Maryland and Rome groups during the SN1987A[niIIllISl- At this 
time the Rome and Maryland detectors showed signals with energy of about 
100 K correlated with the observation of the neutrino detectors Mont Blanc, 
Kamiokande, Baksan and 1MB. Using the previous cross-section, a IQOK signal 
requires the total conversion into gravitational waves of about one thousand so- 
lar masses. This would be impossible for an original star mass of M ^ 20AIsun- 
However, with the renormalized detection cross section of this work, each lOOK 
signal requires just one tenth of a solar mass (still high but energetically possi- 
ble) yielding important information on collapse mechanisms. Also the rennor- 
malized cross-section yields a sensible astrophysical interpretation of the recent 
result concerning experimental coincidences jl4| between the cryogenic resonant 
detectors Explorer (CERN) and Nautilus (Frascati). 

Wc thankfully acknowledge useful comments, discussions and correspondences 
on the subject matter of this paper with Professors N. Cabibbo, E. Coccia, A. 
Degasperis, V. Ferrari, G. Immirzi, F. Palumbo and S. Weinberg. 
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